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We use a toy model to illustrate how to build effective theories for singular potentials. 
We consider a central attractive potential perturbed by a correction. The power- 
counting rule, an important ingredient of effective theory, is established by seeking the min- 
imum set of short-range counterterms that renormalizc the scattering amplitude. We show 
that leading-order counterterms are needed in all partial waves where the potential over- 
comes the centrifugal barrier, and that the additional counterterms at next-to-leading order 
are the ones expected on the basis of dimensional analysis. 



I. INTRODUCTION 



Singular potentials — those that diverge as l/r", n > 2 at a small distance r — are common in 



atomic, molecular and nuclear physics 



l|]. Contrary to the regular sufficiently attractive 



singular potential does not determine observables uniquely Q]. A proper formulation of the extra 
ingredients can be found in the framework of effective field theory (EFT) js], where the problem is 
cast in the usual language of renormalization. With the techniques of renormalization one retains 
the predictive power of singular potentials at low energies, while short-distance physics is accounted 
for in a minimal, model-independent way. 

The basic idea is that in order to obtain the potential from some underlying theory we make an 
arbitrary decomposition between short- and long-distance physics through a short-distance cutoff 
i? ~ 1 /A, above which the potential has a form determined by the exchange of light particles — such 
as two-photon exchange in the case of van der Waals forces, or pion exchange in nuclear physics. 
When the long-range potential is regular, or singular and repulsive, the quantum-mechanical dy- 
namics is insensitive to R. When the potential is singular and attractive, on the other hand, 
observables would depend sensitively on R if we ignored short-distance physics. Dependence on 



the arbitrary scale R can be removed [sl by adjusting parameters of the potential at short dist 



ances 
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as functions of R ( "running constants" ) to guarantee that low-energy observables are reproduced 
independently of R. Such short-range "counterterms" are thus used to mimic the effects of short- 
distance physics. 

The case of the S wave in an attractive singular central potential l/r"' was considered in Ref. 
B, and the particular cases n = 2 and n = 4 were studied in more detail in Refs. Q, B and B, 
respectively. It was shown that a single counterterm — for example, in the form of the depth of 
a square-well — is sufficient to ensure an approximate R independence of low-energy observables, 
once one observable is used to determine the running of the counterterm. An interesting property 
of this running, first noticed in the three-body system [t^, is an oscillatory behavior characteristic 
of limit cycles. For example, when n = 2 the short-distance force has a fixed period as a function 
of Ini?. Limit cycles are now a subject of renewed interest jsj]. 

A consistent EFT must in addition provide systematic improvement over this simple picture. If 
we consider particles of typical momentum Q, we would like to be able to calculate observables in 
an expansion in powers of QRgh = Q/^hi^ where Rsh = is the range of short-range physics. 

Successive terms in this expansion are referred to as leading order (LO), next-to-leading order 
(NLO), and so on. Determining at which orders interactions contribute is called "power counting". 

We are interested here in a two-body system in a regime of momenta where LO comes from a 
dominant long-range potential that is attractive and singular, together with the required countert- 
erms. The latter include one S'-wave counterterm B. However, the long-range potential contributes 
also to other partial waves, and it has been argued on general terms 3] that counterterms are re- 
quired in those waves as well. In the particular case of nucleons, one-pion exchange has an 
tensor force in spin-triplet channels; it has been shown explicitly that counterterms are necessary 
in all waves where the tensor force is attractive and iterated to all orders 
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11 



12]. (For a 



different point of view, see Ref. [131 •) This poses a potential problem because predictive power 
seems to be lost in systems with more than two particles, where all two-particle waves contribute. 
In Ref. |ll| a solution was proposed where, thanks to the centrifugal barrier, perturbation theory 
in the long-range potential is employed in waves of sufficiently large angular momentum. We would 
like to better understand the role of angular momentum, and establish the number of short-range 
counterterms needed in LO to make sense of a singular potential. 

In general a dominant long-range potential suffers corrections that increase as the distance 
decreases. These corrections can arise from additional, smaller couplings to the light exchanged 
particles — e.g. magnetic photon couplings of the Pauli type. They can also be generated by the 
simultaneous exchange of several light particles — such as two-pion exchange. These effects lead to 
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potentials that fall faster at large distances than the LO potential. For example, an LO 1/r" poten- 
tial might have an NLO correction iig^/r""*"^, as is the case for two-pion exchange between nucleons, 
which (neglecting nucleon excitations) goes Q, [l5,[l^ as R'^j^/r^ with Rgh = (47r/^)~^ ~ 0.2 fm 
the characteristic QCD scale. These potentials, which are subleading at distances r^Rgh, are 
more singular than the LO potential and overcome it for r < Rgh- The enhanced singularity ought 
to demand new counterterms. The desired expansion in QRsh requires that NLO counterterms, 
which represent physics at r < Rshi balance the NLO long-range potential — that is, the full NLO 
should be such that changes in observables are small. If that is the case, we expect [lo| that one 
can treat the NLO potential in perturbation theory. One would also expect ll| the number of 
counterterms to be given by naive dimensional analysis — for example, that an NLO correction 
i^^^/r""*"^ needs counterterms with two more derivatives than the counterterms required by the LO 
1/r". In the nuclear case, existing calculations have followed [3] instead the original suggestion 



14 1 that subleading corrections in the potential be iterated to all orders. In fact, it has been argued 



171 ] that this requires fewer counterterms than doing perturbation theory on the corrections. Nev- 



ertheless, treating small corrections in perturbation theory is conceptually simpler, as the running 
of the LO counterterms is not completely modified by the higher singularity of the NLO potential. 
Understanding the strengths and limitations of a perturbative approach would at the very least 
help delineate the scope of a resummation of NLO corrections. 

In this paper we address the issue of power counting for singular potentials, in particular the 
role of centrifugal forces in LO and the perturbative renormalizability of NLO interactions. We 
examine this issue in a simple toy model, where the LO and NLO long-range potentials are taken as 
(central) — 1/r^ and ifcl/r^, respectively. Some of our arguments are similar to those employed in 



related two 



18l | and three [19l . |2C1 . |21| -body contexts. We show that LO counterterms are required 



in all partial waves up to a critical value, and that the number of NLO counterterms is just what 
is expected on the basis of dimensional analysis. Most of our conclusions can be extended to other 
attractive singular potentials. 

This paper is organized as follows. In Sect. [Uwe discuss our EFT framework. In Sect. IIIII a 
new approach for treating the — 1/r^ potential is presented. In Sect. llVl the renormalization of the 
NLO amplitude is discussed and, as a result, NLO counterterms are found. We discuss some of 
the implications of our results to other potentials in Sect. |Vl Finally, we summarize our findings 
in Sect. ED 
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II. FRAMEWORK 

We consider two non-relativistic particles of reduced mass m in the center-of-mass frame, which 
interact through a potential V. Our analysis will mainly be based upon the Lippmann-Schwinger 
equation of the half-off-shell T matrix, 

T{p',p) = V{p',p) + 2mj ^,^^^T{q,P) > (1) 

where p (p') is the initial (final) -state relative momentum. The physics of angular momentum is 
most transparent when we use a partial-wave decomposition. Our convention is that a quantity O 
is given in terms of its partial- wave projection Oi by 

oo 

O(p',p) = ^(2/ + l)O/(p',p)PKcos0), (2) 

1=0 

where is the angle between p and p' . The partial- wave version of Eq. ([1]) is 

m 

Tiip',p) = Viip',p)-- dq ^ \ . Viip',q)Ti{q,p). (3) 

Jo q'^ -p^ - le 

Here we have inserted the ultraviolet cutoff A, which is in general needed to obtain a well-defined 
solution. Since A is arbitrary, observables (obtained from the on-shell T matrix) should be inde- 
pendent of A, 

aJ^T(p,p)=0, (4) 

that is, renormalization-group (RG) invariant. It proves convenient to introduce the reduced partial 
amplitude 

tiip',p)^^Ti{p',p), (5) 
vr 

in terms of which Eq. ([3]) becomes 

tl{p',p) = ^Vi{p',p)- / dq^ ^ -^Vi{p,q)ti{q,p). (6) 

vr^ Jo q'^ — p'^ — le vr^ 

We will also use in numerics the K matrix because of its reality. The K matrix satisfies the same 

integral equation as T but with ie replaced by the principal-value prescription. The reduced form 

for the partial- wave-projected K matrix, 

h{p',p)^^Kl{p\p), (7) 



satisfies 



h{p',p) = ^Vl{p',p)-^ [^dq^^^Vi{p',q)kiiq,p). (8) 



vr^ ' Jo ~ 
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The on-shell ki{p,p) is related to /-wave phase shift 5i{p) by 

tan Slip) fn-\ 

We assume a simple central potential, whose dominant long-range component is singular and 
attractive. We write the potential as 

V{p',p) = Vl{p',p) + Vsip',p) (10) 

in terms of a long-range component Vl and a short-range component Vs- The short-range compo- 
nent Vs{r) is a series of derivatives of delta functions, so Vs{p',p) is a power series in p^, p'"^, and 
p ■ p', starting with a constant. For definiteness, we take the long-range component to be 

VL{r) = vi!'\r) + V^'\r). (11) 

Here V^^ (r) is an attractive inverse-square potential, 

with A > a dimensionless parameter. Its Fourier transform is 

^f(p',p) = -^^, (13) 
■m\p' — p\ 

and its partial- wave projection. 



m{2l + 1) 



p/z+i ' ""^^ ^ I pl+l 



where p> = max{p',p} and p< = min{p',p}. In addition, V^^\r) is an inverse-quartic potential, 

with M a mass scale and g another dimensionless parameter. V^^^ is more singular than vj^\ 
which in momentum space is reflected in higher powers of momenta in the numerator. To define 
the Fourier transform we need to limit the integration to distances larger than a coordinate-space 
cutoff R: 



m {21 + 1) pf+i ' 



(14) 



V^]\p',p;R) = Att / drr^ji{p'r)V^^\r)ji{pr) 
JR 



ir'^g j 4 1 / 21-1 p% 



2mM2 1 ttR {21 + 1){21 - 1) p'-i V 21 + 3 pi 



^'0 + T^^T-W^IPT 1 - + O {RpD . (16) 
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The first term is a constant and cannot be separated from an S-wave sliort-range interaction; it 
can tlierefore be absorbed in V5, and we drop it not to clutter notation. In tlie limit i? — > 0, one 
is tlien left with 

yW(„' - I A (i _ ^LlLA] (u) 

^Ll^P^P)- 2mM2(2/ + l)(2/-l)p'-i V 2Z + 3p2y'- ^ 
We take M to be the characteristic scale of the underlying theory, M = M^i, and \g\ ~ A. In 
this case, vj^^ is a correction to vj^^ at large distances r > l/M^i or, equivalently, small momenta 
p> < Mfii- We want the short-range component Vs to be such that an expansion in Q/Mhi holds 
for observables obtained from the T (or K) matrix. Accordingly, we split Vst T and K as in Eq. 
(jlip with the superscript ^'^^ {^^^) denoting LO (NLO). RG invariance is exact only if all orders are 
considered. Once a truncation to a finite order is made, Eq. can only be satisfied up to terms 
that vanish as A ^ 00 and can be absorbed in higher-order counterterms. In the rest of the paper 
we address the question of which terms should be included in the short-range potential at each 
order to ensure a perturbative expansion of observables consistent with RG invariance. 

III. l/r2 AS LO LONG-RANGE POTENTIAL 

We first tackle the problem in LO, that is, we take V^^^ = 0. The 5 -wave renormalization of 
an attractive inverse-square potential, — l/r^, has been dealt with in coordinate- and momentum- 
space m Refs. 0,Q and Q, respectively. Here we present a new ^jjgroach in momentum space, 
for any partial wave /. Our / = results reproduce known results 
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A. Singularity of 

To see the origin of the peculiarities of a singular potential, we start by taking V^^^ = 0. In this 



case we can write the integral Eq. ([U]) in a simplified form 
ti[px,p) -- 



21 + 1 





e{x-i) - [ 




Jo 



yl + 2 

e{l-x) - I dy^ -ti{py,p) 

— 1 — le 



(18) 



In this form, we see that the only dimensionful parameter is A, and scale invariance is evident in 
the limit A ^ 00 (if it existed). 

The validity of a perturbative expansion in powers of A can be estimated by comparing the first 



FIG. 1: LO amplitude T^°^ as an iteration of the LO potential V'-°\ 



two terms, 

A 



and 



(19) 

21 + 1 ^ ^ 



+ h (20) 



{21 + 1)2 \2l + l 2 

in the expansion of the on-shell ti{p,p) in Eq. ([TS]) . We see that perturbation theory in A works 
well when I is high enough. Conversely, if 

Eq. (I3|) has to be solved for non-perturbatively. The LO amplitude T^^^ is in this case obtained 
from an exact solution of Eq. ([3|) with the LO potential, as shown diagrammatically in Fig. [TJ 

In order to study the non-perturbative regime, we note that the inverse-square potential, Eq. 
(fTSll . has the interesting property 

- V'Vf (p',p) = -±-i2nf6Hp' - p) , (22) 

where V = Using this property we can convert the Lippmann-Schwinger equation, Eq. ([T]), 
into a differential equation, 

V''r(p',p) + \ . Tip',p) = ^5\p'-p) , (23) 

or, decomposed onto partial waves, 

Ax2 



92 

x-Q^ {xtiipx,p)) + 



ti{px,p) = X5{l-x). (24) 



Treating p as a parameter, let us consider the solution of Eq. (p^ . If ti{0,p) is assumed finite, 
which can be inferred from the integral equation, then ti{px,p) for x < 1 is determined up to a 
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coefficient that is a function of p, 

ti{px,p) = Mi{p/A) 2F1 + ^ - ^^«) ' ^ + ^ + '"i^ ' ^ + ^;^^) ' (^^) 

where 2-F'i(a, /3, 7; z) is a hypergeometric function, and 



- (/ + 1/2)2 . (26) 



The pre-factor A// must be a function of p/A on dimensional grounds. It is necessary for calculating 
the on-shell ti{p,p). When x > 1, the solution is the linear combination of two generic solutions, 

1 /, 1 . \ 3 , 2 



ti{px,p) = Mp/A)x' 2F1 1^1 + --mj ,-\^l + - + mj,- + l-x 

with coefficients Ai and Bi that can also be functions of p/A. Since Eq. (j24p is inhomogeneous 
only at X = 1, there is no way to obtain A// unless we determine the ratio of Ai to Bi and then 
match to Eq. (125p . This matching brings cutoff dependence to A//. If the x > 1 solution were cutoff 
independent (up to 0{p'^ /A^) corrections), one could expect Mi and thus ti{p,p) to be constant, 
independent of p — a consequence of scale invar iance. 

The X > 1 solution simplifies in the asymptotic region. For x 3> 1, Eq. (f24|) becomes 

x^{xti{px,p)) + [A - /(/ + 1)] ti{px,p) = O (^^^ . (28) 

This equation has two generic solutions, 

ti{px,p) (xxi-^^'"'") 1 + 



(29) 



If the potential is repulsive, i.e. A < 0, ui is imaginary; one of the two solutions in Eq. (j29p has 
a positive power of x, makes the second integral in Eq. (jlSp divergent when A 00, and has to 
be discarded. When the potential is attractive but not very strong, namely < A < (^ + 1/2)^, vi 
is still imaginary; neither solution produces a divergence in Eq. (jlSp but the one with the bigger 
power must dominate over the other at high x. Hence in these two cases we are left with only one 
solution for x ^ 1, up to a coefficient dependent on p. With the ultraviolet behavior of the off-shell 
ti{p\p) decided one can in principle match to Eq. (j25p . determining thus the on-shell amplitude 
ti{p,p)- 

However, if the potential is attractive and sufficiently strong to overcome the centrifugal barrier, 
that is, A > (/ + 1/2)2, ^j^g^ jg ^gg^^^ pgj. both solutions the second integral in Eq. is finite. 
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but oscillates with A. The two solutions have the same magnitude but different phases. In fact, 
the asymptotic expansion of Eq. ([27j) gives 



ti{px,p) =MI{p/A)x 2 



cos (uiln^ + 61) + ^ cos (ui In^ + 9i + Pi) + ' ^ 



(30) 



A ' 7 ' \ ' A ' ' '-v ' - 

where 9i is an ^-dependent constant and f3i = arg(— 1 + ii'i). Here we inserted the ultraviolet 
cutoff in the cos because the asymptotic dependence should be on p' /A. (We absorbed a factor of 
(p/A)-i/2 in AA'.) The Ion g-range potential, by itself, does not determine the phase 9i + vi Inp/A 
of the solution. In fact, as we change the arbitrary cutoff A, the phase changes. Equation (f30]l is 
the asymptotic form of the solution for p' > p, which is matched to Eq. (I25p to determine ti{p,p). 
Thus, as A changes, so does the on-shell ti{p,p) and the observables obtained from it. 

It is worth noting that the singularity of — depends on the angular momentum I. Equation 
(j26l) implies that for any given A there is a critical Ic above which is no longer real. Therefore 
the singularity exists for 

/ < = \/A- 1/2. (31) 

In each of these waves, the attractive singular potential overcomes the centrifugal barrier, and 
observables are dependent on the arbitrary cutoff. 

We can illustrate these facts with explicit numerical calculations. Since the S wave has already 
been studied in detail elsewhere we focus on the P wave, / = 1. We choose A = 4.25, which is 
strong enough for the singularity to be present in the P wave, that is, 1 < Ic < Ip- Figure [2] shows 
the numerical solution ki{px,p) of Eq. ([8]) for p/m = 0.1, at various values of the cutoff in units of 
the reduced mass, A/m. We see at x ^> 1 the asymptotic oscillations of Eq. ([30l) : both the phase 
and the amplitude depend on A. The A dependence propagates to smaller x and results in very 
different on-shell values ki{p,p). 

Since the choice of A is arbitrary, the LO on-shell k!f\p-,p) should be independent of A when 
A is sufficiently large. But as shown ki{p,p) with V^^^ = oscillates with varying A and the limit 
A ^ cxD is not well defined. Therefore this ki{p,p) and the associated observables are meaningless. 
RG invariance is not being respected in channels where the attractive singular potential is treated 
non-perturbatively and no counterterms are provided. 



B. Renormalization of 

The dependence on A indicates that the T matrix is sensitive to short-range physics through 
virtual states, that is, the loops resummed in the Lippmann-Schwinger equation, Fig. [TJ However, 
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FIG. 2: P-wave reduced K matrix ki{xp,p) as function of x for p/m = 0.1 and various values of A/m, when 
= 0. We have taken A = 4.25. 



the high- momentum part of the loops cannot be distinguished from contact interactions. The cutoff 
simply represents an arbitrary division of where short-range physics is placed. A model-independent 
accounting of this physics requires the introduction of contact interactions with parameters that 
depend on A in such a way that observables are independent of A. Thus, with the technique of 
renormalization one can get a cutoff- independent amplitude. One has to renormalize each singular 
partial wave separately. We shall show here that one counterterm in each singular partial wave 
I < Ic IS necessary and sufficient to renormalize the corresponding partial-wave amplitude. This 
extends the results of Refs. to waves beyond S. 

For the I wave, the contact interaction with fewest derivatives in coordinate space can be written 
in a partial-wave projection as 

.2 r<(o)i 



(32) 



where cj^^^ is a parameter. Equation 



4°-' {px,p) 



becomes 



21 + 1 





e{x-i) - [ 




Jo 



dy 



e{l-x) 



A/p 



dy ■ 



.,1-1 



^ tf\py,p) 



X y 

A/p 



dy ■ 



y 



1-ie 

1+2 







1 



tf\py,p) 



tf\py,p) 



(33) 



instead of Eq. 
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Renormalizing the /-wave amplitude means that the A dependence of C^^^ is such as to make 
the on-shell amphtude tf^\p,p; A^C^'^^A)) independent of A in the large-A hmit: 

aA,<»)(,,,;A,C<»'(A))=0 (£;<!»)) . (34) 

For — 1/r^, we will show that the half-off-shell tf'\px,p] A, C;(A)) is also RG invariant in the large-A 
limit, 

A^tr^{p.,p-A,cr\A))=o(^^,tr^) . (35) 

To justify this claim we take it as an ansatz, evaluate the RG variation of the off-shell tf^\px,p), 
and see if we can make it RG invariant by controlling C^^^ with varying A. To this end, we take the 
total derivative with respect to A of both sides of Eq. (p3]) . assuming that t^^\px,p; A, Ci'^\A)) is 
RG invariant: 

^(f )'{(A-CrA-«)4^A,p) + (pA)W;?fl[l-/,,Kp,A)l} =o(f;r) , (36) 



21 

where we defined 

rA/p 2+2n+l 

InAP, A) = / dy . (37) 

Setting x = A/pin Eq. (j35]) . 

(A,p) = (a - c(°)a-+i) (I)'"^ [1 - IoAp,A)] ■ (38) 

Eliminating the expression in the square brackets in Eqs. ([36]) and (I38p . 



A^ A-C/°)A2W 



21 



(a - C;°)A^'+f - A^('-^^)^] = O . (39) 



Now, the left-hand side of Eq. (j39p can be made of 0{p'^t^^^ / A^) for all x by properly choosing 
C^^'^^(A), which is indeed consistent with the above ansatz that the half-off-shell tf^\px,p) is RG 
invariant, up to 0(p^/A^). This will be so if Ci^\A) satisfies the RG equation 

It is straightforward to solve for Ci^\A), up to a boundary condition: 

^(0),.^ A 2; + l-2;y;tanhln(A/A,;)] 

' ^' A2'+i2/ + l + 2i/;tan[z.,ln(A/A,0] ' ^ ^ 
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FIG. 3: LO P-wave reduced K matrix k'(^\xp,p) as function of x for p/m = 0.1, at various values of the 
cutoff A in units of the reduced mass m, with V^^ introduced to absorb the cutoff dependence. We have 
taken A = 4.25 and A*i/m ~ 0.2. 



mm. 



where A*; is a dimensionful parameter defined in such a way that the reduced coupling 
h^^i^^ cf^\K^i) = —A. It is determined by fitting to the measured I partial-wave amplitude. The 
log-periodic behavior of cf'^ is the so-called limit cycle 

Even with a non-zero {p' ,p) the integral equation ([6]) can be converted into Eq. (j24p . The 
same argument goes through with A^,; substituted for A, dependence on which was eliminated, so 
that the RG-invariant half-off-shell tf'\p' ,p) in the p' ^ p limit is given by 

xp 



tf\px,p)=N'i{p/Ki)x-i 



cos I ui In 



A 



(42) 



instead of Eq. (I30p . The phase is now fixed by A*^, and observables are (nearly) cutoff independent. 

To illustrate this we return to the P-wave example considered at the end of the previous subsec- 
tion. To remove the cutoff dependence observed in Fig. [51 we solve for the K matrix again but now 
with an additional short-range interaction = ir^C^^p'p/Sm put in and allowed to change with 
A in such a way that the on-shell kf^\p,p) is independent of A. Figure [3] shows that not only the 
on-shell k^\p,p) but also the half-off-shell k^i\p',p) are now independent of A, in agreement with 
the previous argument. In this calculation we chose A^,i/m = 0.2, but the results are qualitatively 
the same for any A*i. 

The net effect of renormalization is to replace A by the physical quantities A*; , which parametrize 
short-range interactions and are directly related to data. The appearance of dimensionful param- 
eters (through renormalization with dimensionless parameters and A) is an example of 
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dimensional transmutation. These parameters introduce mass scales in the problem and break 
the scale invariance of the long-range potential — an example of an anomaly. (For an extensive 



discussion of this anomaly, see Ref. [2^.) Because of the log-periodic behavior, however, a discrete 
scale invariance remains, which has striking implications to observables Q]. 

If only one channel (/ = 0) is singular, the trade off between A and A*o is one-to-one; using 
A*o is formally equivalent to treating A as a fit parameter. However, if more than one channel 
is singular at LO, more than one short-range parameter is present in the solution. Ignoring the 
counterterms and just fitting A enforces a link between short-range parameters, which might or 
might not be correct. Regardless of whether this assumption is correct for a given problem, it is 
nothing but a model for the short-range physics, for it is a dynamical assumption that goes beyond 
the symmetry content of the theory. A model-independent treatment of short-range interactions 
requires at LO one counterterm per singular channel where perturbation theory does not apply. 

IV. il/r* AS NLO LONG-DISTANCE INTERACTION 

We now turn to the effects of singular perturbations on the LO singular potential and its 
counterterms. We assume that the NLO long-range potential is given by Eq. (|17p . and ask which 
additional counterterms, if any, have to be supplied at NLO to make the result RG invariant. With 
our choice of parameters, the long-range potential, Eq. (jl7p . is a correction to the LO long-range 
potential, Eq. (I14p . so the full NLO should have a small effect on observables. Accordingly, we 
treat the NLO potential V^^^ in perturbation theory, that is, in first-order distorted-wave Born 
approximation. The NLO amplitude T^^) has one insertion of the NLO potential, see Fig. [H that 
is, 

fA 



^ Jo q — p — te TT 

/"A /"A ^/ /my-,' 



tf\q,p). (43) 



— p"^ — ie ' 

where we have used the symmetry of Vi under exchange of its arguments. Two insertions of V^^'^ 
(second-order perturbation theory) come at next order, where further contributions to the long- 
and short-range potentials might exist. 

For simplicity, we focus on the 5 wave, / = 0, in the following. Generalization to higher waves 
and higher orders is straightforward. 
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FIG. 4: NLO amplitude T^^^ in first-order distorted- wave Born approximation in the NLO potential F*^^^. 

A. Additional singularity 



In principle new cutoff dependence, possibly divergent in the A ^ cxd limit, will arise from the 
loops in Fig. 01 Omitting for the moment any additional contact interactions in NLO, that is, 
taking V^'^ = 0, let us consider the A dependence of the on-shell amplitude. Taking the total 
derivative of Eq. with respect to A and using the fact that the off-shell t^Q^{jpx,p) is RG 

invariant up to suppressed terms (see Eq. ([35|) ). we find 



A^to(p,p) 



(l {1 - ^o,o(p, A) + 1^ [1 - A,o(p, A)] \ tf{Kp) 



+0 



AO) 



VM2A2 

where the integrals In.i{p,-^) were defined in Eq. ()37p . Using Eqs. ()38p and ()42p we can write 
[1 - /o,o(p,A)]4°)(A,p) = AAo2(p/A.o) k + (A\+o(^^ 



(44) 



(45) 



where Go and Gi are two p-independent, oscillating functions of A/A*o- A similar form can be 
obtained for the other term, 



p2_ 

A2 



[l-/i,o(p,A)]4°)(A,p)=AAo2(p/A,o 



■'2/ 



A 



+ A^""^ I a:^ 



+0 



(46) 



A/p /-A'/p 
' 2,(0), 



in terms of two other oscillating functions Hq and Hi that do not depend on p. This can be seen 
from the cutoff dependence of Iifi{p,A) in the large-A limit. We can write 

dyy^t^o\py,p) , (47) 

where A' is a scale above which the asymptotic expansion 

of 4°\ Eq. g2]), is valid. Since the A 
dependence of the second integral is at most oc p^ j (see Eq. ([35]) ). Hq and }i\ are given by the 
first integral. 



/i,o(p,A) = /o,o(p,A) + / dyy't'^\py,p) + 

lA'/p JO 
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Thus, Eq. 



can be written as 



A^to(p,p) 



A. 



*o 



+ p 



A 



A: 



*o 



A 



A, 



*o 



A 



Gi{^] +Go[ — ] +-Hi[ — ] +-Ho[ — 



A: 



*o 



1 



A 



A 



*o 



We find that there are two types of A-dependent terms (modulated by AAq^) that do not vanish 
in the large-A limit: an energy-independent term whose oscillatory behavior gets enhanced by 
an arbitrary factor A^, and a term that introduces further cutoff dependence proportional to 
the energy. This stronger cutoff dependence is just the momentum-space reflection of the higher 
singularity of Results become sensitive to the physics at the smaller distances where vj^^^ 

overcomes Vj^\ To account for this physics in a model-independent way, new counterterms are 



needed. 



B. NLO renormalization 

Since is more singular than vj^^ by two powers of momenta (c.f. Eqs. (I14p and ()17p ). we 
expect, on the basis of dimensional analysis, that new counterterms with up to two extra derivatives 
or powers of A will be required. Indeed, the two types of A dependence in Eq. (I48p suggest that we 
need two new counterterms, one being possibly just a correction to In that case, the running 
of Co is changed at NLO. For clarity, we split Cq into two pieces, Co = Cq^^ + Cq^\ where Cq'^^(A) 
has the LO running given in Eq. gl]) and C^^\A) has an NLO running to be determined. However, 
this counterterm cannot be expected to eliminate the energy-dependent term. That requires a new 
counterterm D^^ , which represents the leading energy dependence of the short-range physics, and 
whose running Dq^^ (A) should also be determined from the requirement of RG invariance of the T 
matrix. These arguments suggest that the NLO short-range potential is 



m L 



(49) 



Including these terms and using Eq. ([38]) . we get, instead of Eqs. (jH]) and (08]), 



A2 [1-/o,o(p,A)] 



+ p2 [i-/i^o(p,A)] 




+0 



( P 



VM2A2 



-N'q^{p/A^q\ 



A^ 

M2 



R{A) + ^S{A) + o(^ 



(50) 
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where 




and 




9 |0^(i)a 2A^ ^^o'^ ' 
.(1) 



It is clear that C^p and generate terms with the same types of A dependence as in Eq. 
Therefore, by controhing Cq^^ and D^q \ we can make both -R(A) = and S{A) = 0, and the NLO 
amphtude t^Q\p,p) becomes RG invariant up to 0{p^ /M"^!^). 

Because the RG equations of the NLO counterterms Cq^''(A) and D^q\a) are not particularly 
illuminating, we turn to numerical experiments. We take A = 2, g( = 1 and M = 0.5m. To test RG 
invariance, we define the fractional NLO correction 

4\p,p;A) 



X{p,A) 



kf\p,p;A) 



(53) 



We first show that, in agreement with the previous argument, renormalization cannot be done 
with Cg^'* alone. We take as "datum" A;o(0.1?ti, 0.1m) = —1.05. In LO, C^\a) is determined so 
as to reproduce this datum. The energy dependence is a prediction of the theory. In NLO, the 
additional terms in the potential will make the theory deviate from the datum, unless Cq'\A) is 
fitted to preserve agreement with the datum. We thus solve the Lippmann-Schwinger equation 
with various cutoffs, and find C^\a) such as to yield the datum. In Fig. [5] the dot-dashed line 
shows the fractional NLO correction X[p,A) as function of A for p = 0.175m. It oscillates as 
A increases and does not show sign of convergence. One concludes that by itself does not 
renormalize the NLO amplitude. 

We now repeat the calculation but including both and Since two parameters are 

involved in the fit, /co(0.1m, 0.1m) = —1.05 and fco(0. 15m, 0.15m) = —0.34 are used as "data". 
The second data point is chosen as a 5% displacement of the value of fcp'^^ (0.15m, 0.15m), to ensure 
that at low momenta NLO represents a small effect on observables. The cutoff dependence of the 
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FIG. 5: Fractional NLO correction A) at p = 0.175m as function of A/to. The dot-dashed line is made 
by setting Dq^' = while the solid line employs both Cq^'' and d'^^^ ■ We have taken A = 2, y = 1, and 
Mjra = 0.5. 

fractional NLO correction at p = 0.175m is shown as the solid line in Fig. [5j The plateau in the 
solid line supports the hypothesis that the counterterms in Eq. (I49p indeed renormalize the NLO 
amplitude. 

With the data chosen above we compute the energy dependence of the amplitude. Figure [6] 
shows the energy dependence of both LO and LO+NLO reduced K matrices. Both sets of results 
are computed with four different cutoffs A/m = 5.5, 6.5, 7.5 and 8.5: there are four data points on 
each spot in the figure. The fact that the points for different cutoffs are indistinguishable indicates 
that the amplitudes are being properly renormalized. It is seen that the NLO correction is small 
for p <^ M and fails for p Af , as it should. 

V. DISCUSSION 

In the previous sections we tackled a number of issues in the rich physics of singular potentials 
using a simple toy model — 1/r^ it 1/r^. Most of our results are independent of this particular 
choice, but some are a consequence of the classical scale invariance of the LO potential. We discuss 
here some of the limitations and generalizations of these results to an LO potential of the form 
-A/2mAr,""V", n > 2 and A = 0(1) > 0, perturbed by a more singular potential. Here Mio 
is a characteristic scale that sets the curvature of the long-range potential. Clearly, at distances 
r ~ 1/Af/o there can be a balance between kinetic repulsion and long-range attraction, so this is 
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FIG. 6: Reduced on-shell K matrix ko{p,p) in the S wave as a function oip/m, with A/to — 5.5, 6.5, 7.5 and 
8.5: leading order (crosses) and up to next-to-leading order (circles). We have taken the same parameters 
as for the solid line in Fig. [3 



a natural size for a bound state or resonance to have. Accordingly, a momentum Q ~ M/o is the 
most interesting resolution to consider. 

In Sect, mil we presented a new method for the study of the Lippmann-Schwinger equation with 
a singular potential, in which we transformed it to a differential equation, Eq. (|23p . We would like 
to point out that this method can be applied to all 1/r^™ (m an integer) potentials, in which case 
we obtain a generalization of Eq. ()23p that involves the mth-order Laplacian ^ . The corresponding 
solutions then have similar properties. An extension of the calculations of this paper to the more 
general case is under investigation [23]. 



We have shown that angular momentum plays an important, double role through the repulsive 
effects of the centrifugal barrier. The emergence of the remarkable phenomena we discussed comes 
from the competition between the singular potential and the centrifugal barrier. When n = 2, 
for a given singular potential strength, there is a critical angular momentum Ic, Eq. ()3ip . above 
which the effective radial potential is no longer attractive, and the particles are prevented from 
probing short-range physics. In these waves, the problem is well defined in LO without a short- 
range counterterm. Conversely, in waves with I < Ic, a counterterm is required in every wave, in 
agreement with a general argument [9]. Such a critical angular momentum is a particularity of 



^ The resultant differential equation resembles a Schrodinger equation, since the operator is just the Laplacian 
in the momentum representation. It is not clear to us if a similar transformation exists for any type of long-range 
potential. 



19 



n = 2. For n > 2, the singular potential will overcome the centrifugal barrier (with / > 0) at 
some distance ri ~ (A//(/ + /Miq. Therefore, the two particles will want to get as close 

as allowed by short-range physics. To make observables minimally sensitive to this short-range 



11 



1^ 



physics, a single counterterm is required in every wave. Only recently was this found 
in the more complicated case of the direction-dependent 1/r'^ potential originating from one-pion 
exchange between nucleons. 

Regardless of the existence of Ic, the centrifugal barrier weakens (for n = 2) or dominates (for 
n > 2) over the singular potential at large distances. The larger the I, the weaker the effective 
radial potential in a partial wave. This leads for n = 2 to a second value for /, Ip in Eq. (j2ip . 
above which the potential can be treated in perturbation theory. When n > 2, there is always a 
spatial region r < r/ where the singular potential is strong, but its size decreases with /. Whether 
the potential can be treated in perturbation theory in a given partial wave depends then crucially 
on the range of momenta Q that we are interested in probing. For a given Q, a sufficiently large / 
exists where the support of the strong potential is effectively a short-range effect. For short-range 
potentials, barring fine-tunings ^ leading to bound states or resonances at threshold, perturbation 
theory holds in higher waves In the nuclear case below the QCD mass scale, it was found that 



Ip 



111. 



It is thus a general feature of singular potentials that a single counterterm, Eq. (|32p . is needed 
in all waves where the potential is sufficiently attractive. For n = 2, the presence of dimensionful 
parameters A^,; breaks scale invariance, but discrete scale invariance remains in the form of limit 
cycle. Our LO results in the S wave are essentially the same as in the three-body system with 
short-range interactions . For n > 2, the dependence of counterterms on the cutoff is slightly 

more complicated, but still oscillatory 3, ^. 

The renormalization of — 1/r" suggests a different power counting than what one would expect 
from naive dimensional analysis (NDA) or naturalness. Based on NDA one would expect that Ci 
in Eq. ([H^]) scales as 

(-) 

This is reasonable if A is outside of the singular region or if the Born approximation is valid. But if 
this is not the case the necessity of renormalization promotes C/ to the same size as the long-range 
potential. The renormalized C/ can be thought of as scaling instead with the resolution Q at which 
the potential is considered valid. 



C-r-Q^. (55) 
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The corresponding short-range interaction, Eq. (|32p . then scales as 1/mQ, which is the same 
scahng as the long-range potential. 

In Sect. llVl we investigated the effects of a perturbation in the form of a singular potential with 
two more powers of momenta than the LO potential: a {Q/Mhi)"^ perturbation. In the nuclear case, 



it has been suggested 



111 and disputed 



171 ] that the additional counterterms are those indicated 
by NDA, in this case, those with two more powers of momenta. We found here that a perturbative 
treatment of the corrections is indeed consistent with this NDA expectation. 

Therefore, once the failure of NDA is corrected at LO, Eq. (j55p . power counting is formulated 
as usual. In the particular example considered here, Eq. (|49p . 



- 4- ' - ^ ■ (56) 

° Ml' ° QMl ^ ^ 

As a consequence, the NLO short-range interaction, Eq. (j49p . scales as gQ /mMl, just as the NLO 
long-range potential. 

We can also use the arguments of Sect. IIVI to determine the NLO counterterms if the NLO 
long-distance potential is another potential than 1/r^. For example, in the case of as 
NLO, whose Fourier transform is ~ Q"^ /niMl, one expects to need 

ci'^+^r(/+p') + 4VV (57) 
as NLO counterterms in the S wave, based on the qualitative arguments of Sect. IIVI In this case 

C^'^-'X, Dl^^r^aQ E^^^^^, (58) 

so that the short-range interaction scales as gQ^ / mM^-. 

This state of affairs is perhaps not surprising. NDA was developed based on experience with 
perturbation theory. It does fail for an attractive singular potential, but only when the potential is 
treated non-perturbatively in LO. Once this unforeseen event is incorporated in the power counting, 
the perturbative treatment of the corrections conforms to NDA, relative to the corrected LO. 



VI. SUMMARY AND OUTLOOK 

We have used — l/r^(LO) it l/r^(NLO) as an example to demonstrate how to build effective 
theories based upon singular potentials. The key point is to understand the power counting of 
contact interactions. The correct power counting scheme should consist of the minimum set of 
contact interactions that renormalizes scattering amplitudes including the long-range potentials. 
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Due to this intrinsic relation between renormalization and power counting it was found that the 
sizes of contact interactions are different from what is expected by naive dimensional analysis. 

A new approach to the renormalization of a — 1/r^ potential was presented. It was shown 
that the singularity of — 1/r^ depends on angular momentum. The region of singularity entangles 
with that of validity of the Born approximation. In each partial wave where the LO potential is 
resummed to all orders, a single counterterm is needed for renormalization. The NLO potential 
can be treated as a (distorted-wave) perturbation, and the minimum set of NLO short-range 
counterterms that are needed to renormalize the NLO amplitude can be determined by estimating 
the superficial cutoff dependence with the asymptotic LO T matrix. Analytical arguments were 
supplemented by numerical evidence. 

It is one of the main conclusions of this paper that the power counting, a key ingredient of any 
effective theory, cannot be decided solely on the basis of naive dimensional analysis in the case of 
(non-perturbative) singular potentials. One has to rely on explicit checks of RG invariance, either 
analytical or numerical, in order to test any proposed power counting scheme. Yet, once LO has 
been understood, perturbative corrections do not violate naive dimensional analysis with respect 
to LO. 

3,Q, 



Besides nuclear forces 



101, 



11 



12 



151, 



161, 



13], there may be other applications of effective 



theories of singular potentials. For example, the 1/r^ potential in two dimensions is relevant for 



251]). while it appears in three 



the interaction of a neutral atom with a charged wire (see, e.g., Ref. 
dimensionswith an angle-dependent coefficient when a charge is in the field of a polar molecule (see, 
e.g. Ref. [26]). Our results could be readily applied to long-range corrections in these systems ^. 
With extensions I23], it could also be applied to the (long-range) electron-atom interaction, which 
is often divided into — 1/r^, 1/r^ and higher terms [1]. In all these cases, one can construct effective 
theories with a well defined power counting that incorporate renormalization-group ideas. 



^ It is well known that a three-body system with short-range pairwise interactions can be mapped [231 — when the 
two-body S-wave scattering length a2 00 — into a two-dimensional — potential problem. Our method can 
thus be adapted to this system as well. The counterterm necessary to renormalize the — problem in LO 
represents a three-body force in the original system, while our NLO analysis is related to the controversy [2^ 
regarding the renormalization of higher-order three-body forces in three-body systems when both 02 and the two- 
body effective range r2 are finite. However, if 02 is finite the mapping is complicated [2l| even in the case r2 = 0. 
An investigation of whether our method can be usefully applied to the three-body case of interest is worthwhile 
but beyond the scope of the present manuscript. 
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